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Abstract 

The dependence on the single-particle states of the pairing matrix 
elements of the Gogny force and of the bare low-momentum nucleon- 
nucleon potential viow-k - designed so as to reproduce the low-energy 
observables avoiding the use of a repulsive core - is studied in the 
semiclassical approximation for the case of a typical finite, superfluid 
nucleus (^^''Sn). It is found that the matrix elements of viow-k fol- 
low closely those of vcogny on a wide range of energy values around 
the Fermi energy ep, those associated with viow-k being less attrac- 
tive. This result explains the fact that around ep the pairing gap 
^Gogny associated with the Gogny interaction (and with a density of 
single-particle levels corresponding to an effective /c-mass ~ 0.7m) 
is a factor of about 2 larger than Aiow-k, being in agreement with 
^eap= 1-4 MeV. The exchange of low- lying collective surface vibra- 
tions among pairs of nucleons moving in time-reversal states gives rise 
to an induced pairing interaction Vind peaked at ep- The interac- 
tion {vio^-k + Vind)Z^ arising from the renormalization of the bare 
nucleon-nucleon potential and of the single-particle motion (cj— mass 



1 



and quasiparticle strength Z^^) due to the particle- vibration couphng 
leads to a value of the pairing gap at the Fermi energy ^^^^i which 
accounts for the experimental value. 

An important question which remains to be studied quantitatively 
is to which extent Ac^gny, which depends on average parameters, and 
Arem which explicitly depends on the parameters describing the (low- 
energy) nuclear structure, display or not a similar isotopic dependence, 
and whether this dependence is born out by the data. 



1 Introduction 

An economic description of pairing correlations in finite nuclei is provided by 
Hartree-Fock-Bogoliubov (HFB) theory ^1 12] making use of phenomenologi- 
cal interactions like e.g. the finite range Gogny force 3J or density dependent 
zero-range forces combined with appropriate energy cut-offs (cf., e.g., jHIHl)- 
Such a description leads to values of the pairing gap which are in overall 
agreement with the experimental findings. Note that in these calculations, 
the density of levels at the Fermi energy piep) is controlled by the so called 
fc— mass (i.e. piep) ~ mk) which, as a rule, is smaller than the bare nucleon 
mass (e.g. ^ 0.7m in the case of the Gogny force). 

On the other hand, a number of studies have shown that the super- 
fluid properties of nuclear systems, ranging from infinite nuclear and neu- 
tron matter to finite atomic nuclei, are strongly influenced by polarization 
phenomena 6^ . In these calculations one starts from a bare nucleon-nucleon 
potential (Argonne, Bonn, Paris, etc.) adding afterwards the renormaliza- 
tion processes. Recently, the pairing gap, the quasiparticle spectrum and the 
collective modes of ^^°Sn have been calculated solving the Dyson-Gor'kov 
equation j7j, in a single-particle space characterized by = 0.7m, allowing 
the nucleons to interact in the ^5*0 channel through a f 14 Argonne N-N poten- 
tial taking into account the variety of renormalization processes (self-energy, 
fragmentation, induced interaction and vertex corrections) arising from the 
coupling of the particles with surface vibrations. While the bare N-N inter- 
action accounts for about half of the pairing gap, overall agreement with the 
experimental findings is achieved by including medium polarization effects. 
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In the present paper we want to shed hght into the physics of these results, 
by studying the magnitude of the different pairing ( J'^ = O"*") matrix elements 
as well as their dependence on the energy of the single-particle states lying in 
the vicinity of the Fermi energy in Section 2, and by discussing the calculation 
of the pairing gaps associated with these matrix elements in Section 3. 

2 Matrix elements 

The matrix element of the induced interaction can be written as jH] 



where L, M and n denote the quantum numbers of the exchanged collective 
vibrations, fLnij) = PLnRo{dU/dr) is the associated radial formfactor is 
the deformation parameter associated with the n— th mode of multipolarity 
L, Rq is the ground-state radius and U is the average potential), hu^n is the 
vibrational energy of the n— th mode, while and ep are the single-particle 
and Fermi energies. Eq is the pairing correlation energy per Cooper pair, 
which is of the order of —A. In practice we have used Eo= - 2 MeV. The 
pre-factor 2 in Eq. (jT)) comes from the two possible time orderings associated 
with the one-phonon exchange. 

In Fig. Qwe show the value of the diagonal matrix elements (z/z/|f j„(^|z/z/) 
of the induced interaction for the nucleus ^^°Sn as a function of the single- 
particle energy. The single-particle levels entering Eq. have been calcu- 
lated making use of a standard Saxon- Woods potential, whose parameters 
are provided in the caption of Fig. ^ In most of our calculations we use this 
potential and an associated effective mass equal to the bare mass. In the 
last section we shall use instead an effective mass m^, = 0.7m, which simulates 
the typical outcome of a Hartree-Fock calculation (with, e.g., the Gogny in- 
teraction). In those cases, the Woods-Saxon parameters are changed in such 
a way that the Fermi energy remains close to the experimental value (V^ 
^ (m/mfc)Vo 0, see e.g. Fig. |SJ. The parameters Pin and huLn have been 
determined by diagonalizing a separable multipole-multipole interaction in 
the quasiparticle-random phase approximation (QRPA) adjusting the cou- 
pling constants so as to reproduce the energies and transition probabilities 
of the lowest-lying states of each spin and parity. Also shown in Fig. 1 are 



LMn 




(1) 



3 



the diagonal matrix elements of the Gogny interaction. Here, and in the rest 
of the paper, we shall employ the Dl parametrization. 

The matrix elements of Vind-, which display an average value equal to -0.15 
MeV, are peaked at the Fermi energy, within an energy range of essentially 5 
MeV around ep- This behaviour reflects the fact that (z/'z/'|wj„rf|z/, i/) is con- 
trolled, through the energy denominator appearing in Eq.([Q), by the energy 
of the low-lying collective vibrations {huiin) < 5 MeV). Outside this energy 
range, the values of (i^'z^'|fmd|^, ^) become very small, less than 50 keV in ab- 
solute value. The matrix elements of vcogny decrease on average smoothly in 
magnitude all the way from the deepest levels up to the continuum threshold, 
displaying a value of about -0.3 MeV ep. 

The induced interaction matrix elements depend on specific properties 
of finite nuclei, namely the single-particle quantum numbers and energies, 
and the energy, the zero-point amplitude and the transition density of the 
vibrational states. To study the global features of the bare, of the induced, 
and of the Gogny interactions as well as to provide a more transparent form 
for their matrix elements, we shall present, in the rest of this section, also 
the matrix elements obtained using the semiclassical Thomas-Fermi (TF) 
approximation. In fact, this approximation averages out shell effects typi- 
cal of finite nuclei and provides the overall energy dependence of the matrix 
elements. Within this context we follow Ref. ^U] where semiclassical expres- 
sions for two-body matrix elements have been presented. To obtain a smooth 
behaviour of the diagonal matrix elements as a function of the continuous 
energy variable we compute the Fourier transform of a two-body inter- 
action f(ri,r2) using plane waves of momenta pi and p2- The average of 
the quantal matrix elements is carried out using the normalized semiclassical 
density matrix. 



Rc being the classical turning point. We assume a constant value for the effec- 
tive mass rriki neglecting its possible dependence on position. The expression 
for the semiclassical diagonal matrix element of f (ri, as a function of the 





In the above relation the semiclassical level density g is 
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single-particle energy, is consequently (TU] 



v{E) = J d rid r2v{fi,f2) x 

d^Pid^P2e^^P'-P''^\ pI pI 

i2nnr g\E) -2^,- 

For an interaction that depends only on the magnitude of the relative co- 
ordinate s = r*! — r2, like the Gogny interaction in the S=0, T=l channel, 
one can integrate first over d^s and then over the momenta d^pi and d^p2, 
obtaining the expression 

v{E) = c{E) J d^R{E - U{R))e{E - U{R))v{kE{R), kE^R)), (5) 

where R indicates the center of mass, v{k, k) is the on-shell pairing matrix 
element of the interaction v{s) in the ^5*0 channel and in a plane wave basis. 



evaluated at the local momentum k = kE^R) = ^y2mk{E — U{R)), and 
c{E) denotes the quantity 

In the next section we shall calculate the pairing gap associated with the 
different interactions, and therefore we shall need the matrix elements be- 
tween all pairs of particles moving in time reversal. For this purpose we shall 
employ a formula analogous to that given in Eq. (0), namely 

v{E,E') = ^c{E)c{E') J d^R^{E - U{R))^{E' - U{R))x 

X e{E - U{R))e{E' - U{R))v{kE{R), kE'{R)). (7) 

In this relation E and E' indicate the energies of the single-particle states of 
each pair. 

We now turn to the case of the induced interaction, associated with the 
matrix elements of Eq. (1). In this case the matrix elements depend sepa- 
rately on fi and r2. We then integrate over the momenta first, obtaining 

VindiE) =2}^ C{E) I I 7- : 

(8) 
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where jo is the Bessel function and k = y -p^(i? — U{R)). The index n labels 
the different phonons of multipolarity L (the phonons included in the present 
calculation are specified at the end of this Section). Note that one can obtain 
numerical results from Eq. (jH)) by avoiding the six-dimensional integration. 
This is achieved by performing a multipole expansion of the quantity 

A'TT 

jl{ks){E - Um = Y: ^j—YFiin, r^, E)YUri)YiUf2), 

Im 

which leads to 

^ Eo-2\E - eFl-hiULn 

In Ref. ^] the reliability of the TF approximation in reproducing the 
quantum mechanical (QM) matrix elements has been checked in the par- 
ticular case of the pairing matrix elements of a 5-interaction acting among 
particles moving in a harmonic oscillator (HO) potential. Note that the TF 
approach implies an averaging over the quantum numbers associated with 
the different states belonging to the major shells and displaying an energy 
En = {N + ^)hu!. The TF matrix elements have been compared with the 
averaged QM matrix elements (taking degeneracies into account). Overall 
agreement was observed (cf. Table II of Ref. [TU]). 

A similar comparison for the pairing matrix elements of the Gogny force 
in the ^5*0 channel is reported in Fig. |21 It is seen that the TF approach pro- 
vides a good approximation to the average quantum matrix elements. How- 
ever, the QM diagonal matrix elements, in particular those associated with 
s-states, are systematically larger (in absolute value) than the non-diagonal 
ones. One should be aware of the fact that the semiclassical density matrix 
can be considered to be the analogue to the one obtained from a Strutinsky 
smoothing. Therefore, Eq. (0), implicitly, represents a function of about Ihuj 
width. Thus, the matrix element in Eq. (0} contains at the quantum level 
cross terms corresponding to at least one major shell. Taking for example the 
2s-ld shell, we have to weigh the two diagonal matrix elements with the 2s 
and Id wavefunctions with a factor 1 and 25, respectively. The non-diagonal 
2s-ld element obtains the weight factor 10. Performing the arithmetic aver- 
age yields the crosses of Fig. |21 In this way, we see that the TF expression 

reproduces very well the quantal average. The more realistic case of the 
Gogny matrix elements in the Woods-Saxon potential (including spin-orbit) 
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is illustrated in Fig. |21 The general pattern is similar to that shown in the 
previous figure. 

Having assessed its validity, we now employ the semiclassical approxima- 
tion to compare the matrix elements of the induced interaction, the effec- 
tive Gogny interaction and a bare nucleon-nucleon interaction. Concerning 
the latter, the matrix elements of bare nucleon-nucleon interactions which 
contain a repulsive core display a qualitatively different momentum depen- 
dence than that displayed by the Gogny interaction. We shall instead con- 
sider here the viow-k interaction, which has been devised to reproduce the 
low-momentum properties of the nucleon-nucleon force, including the exper- 
imental phase shifts, without the introduction of a repulsive core ^Tj. In 
particular, we shall employ a parametrization of viow-k devised to reproduce 
the properties of V14 Argonne potential. This parametrization requires a mo- 
mentum cutoff kcut equal to 2.1 fm~^. The matrix elements of viow-k in a 
plane wave basis have been calculated in Ref. 12 and are reported in Fig. |3] 
It is seen that, while Vcogny is less attractive than viow-k for small momenta, 
it becomes more attractive for k larger than ^1 fm~^. 

The semiclassical matrix elements of the Gogny Dl interaction (Eq. (j3j)) 
and of the viow-k interactions calculated as a function of the single-particle 
energy for the nucleus ^^''Sn, are displayed in Fig. |S1 We remark that, close 
to the Fermi energy, the matrix elements of the Gogny interaction calculated 
with the DIS parametrization are about 10 keV less attractive. It is seen 
that the general trend is the same as that obtained in infinite matter, with 
the difference that the matrix elements of the two interactions look somewhat 
more similar to one another. Insight into this difference can be obtained with 
the help of Eq.(jS)). This relation is essentially an average of v{k) calculated 
at the various local momenta in the nucleus. From Fig. |3]we see that in the 
interior of the nucleus (at larger local momenta) Vcogny is more attractive 
than viow-k, but on the surface it tends to be less attractive than viow-k- 

The semiclassical matrix elements of the induced interaction are also 
shown in Fig. |31 These matrix elements are strongly peaked around the 
Fermi energy and reproduce rather well the average of the quantal results al- 
ready shown in Fig. ^ In Fig. [Ulwe have separated the contributions arising 
from the exchange of phonons of multipolarities L=2, 3 and 4. We include all 
the phonons with energy below 30 MeV in the calculation. The importance 
of the contributions with L=2 and 3 can be understood in terms of the high 
collectivity and low energy of the associated collective modes. The energies 
and deformation parameters of the low- lying modes are huj2 = 1.17 MeV , 
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hujs = 2.42 MeV, ku^ = 2.47 MeV and /^a = 0.12, (3^ = 0.15, = 0.07. 

We also note that at the Fermi energy the matrix elements of Vind are 
rather large, of the order of those associated with the bare nucleon-nucleon 
potential viow-k (in this connection, cf. the estimate in ref. [T3]) . 

3 Pairing gaps 

In this section, we shall discuss the results of calculations of the pairing gap 
which employ the matrix elements discussed above. In the nucleus under 
study, namely ^^°Sn, we obtain the state-dependent pairing gaps by solv- 
ing the generalized BCS equations [EQ: this means that we solve the HFB 
equations by treating the pairing sector self-consistently, while the mean 
field is described in terms of the Woods-Saxon potential already used in the 
calculation of the matrix elements. In this way, we take into account scat- 
tering processes between nucleons lying on orbits having different number 
of nodes/energies associated with non-diagonal matrix elements of the 
type v{Ei, E2; E[, E2), with different energies Ei and E2 (or E[ and E2). 
This is at variance with the usual BCS method in which the pairing gap 
receives contributions only from scattering processes between time-reversal 
states. We shall refer in the following either to quantal or to semiclassical 
calculations; in the latter case, this is done by solving the same equations 
as in the quantal case, but replacing the QM matrix elements with the TP 
ones, using Eq. ((Tj) and the prescription 

E={Ei + E2)/2, E' = {E[ + E'2)/2. (10) 

The positive energy states are obtained by setting the system in a spheri- 
cal box. We have checked that convergence is achieved using i?{,ox= 12 fm and 
including all the states from the bottom of the potential up to the positive 
energy Ec«t=50 MeV. 

In infinite neutron matter, the pairing gaps calculated with the interac- 
tions viow-k and VGogny are very similar up to Fermi momenta of the order of 
0.6 fm~^, becoming increasingly different from each other at higher momenta, 
as can be seen from Fig. [7| (cf. Ref. [HIE], where the DIS parametrization 
of the Gogny interaction has been used). The pairing gap obtained with 
viow-k goes to zero, around saturation, as a function of kp much faster than 
that associated with Vcogny] it can also be seen that the pairing gap obtained 
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with viow-k reproduces quite accurately the result found allowing particles to 
interact through the Argonne V14 bare N-N potential [T^ . 

In the case of ^^°Sn, the pairing gap with the V14 Argonne interaction 
and with = m, was calculated previously in ref. obtaining a value 
of about 2.2 MeV close to the Fermi energy. We have verified that essen- 
tially the same value is obtained with viow-k, in keeping with the infinite 
matter case (cf. Fig. [7j). In the following, however, we shall adopt the 
value TUk = 0.7m, which is the typical value associated with Hartree-Fock 
calculations performed with Gogny or Skyrme forces. In Fig. |H1 we compare 
the state-dependent pairing gap obtained with the Gogny force and with the 
viow-k interaction. We also compare the pairing gaps obtained inserting the 
semiclassical pairing matrix elements in the HFB equations. We remark that 
the semiclassical approach takes out the scatter from shell effects which is a 
desirable feature when one wishes to reveal generic trends. 

The Gogny interaction leads to a pairing gap of 1.4 MeV close to the 
Fermi energy, in good agreement with the experimental value deduced from 
the odd-even mass difference. This value is about two times larger than that 
obtained with viow-k- Note also that the pairing gap obtained with viow-k 
reproduces the results obtained with the bare W14 interaction |7j . We remark 
that using the DIS parametrization of the Gogny force, instead of Dl, one 
obtains a gap at the Fermi energy of about 1.2 MeV. 

In Fig. ini we show the state-dependent pairing gap obtained with the 
pairing matrix elements of Vtotai = Viow-k + Vind, resulting from the sum of the 
matrix elements of the bare interaction viow-k and of the induced interaction 
Vind- The average value of the resulting pairing gap ^totai at the Fermi 
energy is about 2 MeV. This value is about 25% higher than that obtained 
in Ref. [7|, where we solved the Dyson-Gor'kov equation, taking properly into 
account all polarization effects (induced interaction, self-energy and vertex 
corrections) |17|. 

This is consistent with the fact that according to Nuclear Field Theory 
(NFT) PHI; if one considers the effects of the exchange of phonons between 
pairs of nucleons, one has to consider at the same time processes where 
the phonon is absorbed by the same nucleon which has virtually excited 
it. Such processes lead to self-energy (w-mass and single-particle splitting) 
as well as to vertex correction phenomena. Of these processes and for the 
nucleus under discussion (^^°Sn), u;— mass effects P are the most important. 
Taking such effects into account is equivalent to use a residual interaction 
Vren = VfotaiZoj, whcrc = {m^ / m)~^ is the quasiparticle strength jU] at 
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the Fermi energy (cf. Appendix A). Consequently Aj.e„ = Z^Atot- Because 
Vind is proportional to Pi/ DENl, that is to the square of the deformation 
parameters associated with the low-lying collective vibration and inversely 
proportional to the energy denominator DENl appearing in Eq.(^, a typical 
error of 20% on Pl and DENl implies a 60% error in Vmd- Making use of 
this fact and of the results of Figs. [T] and IHl one obtains, as discussed in 
Appendix B, Aren ~ Z^\ot ~ 1-35 ± 0.15 MeV. It is seen that the 60% 
error in v^d has been reduced to a 10% in Aind- This is because the stronger 
the particle- vibration coupling is, the larger Atot but the smaller the amount 
of single-particle content of levels around ep (and thus the smaller Z^^), and 
viceversa. The self consistency between collectivity of the modes, strength of 
the particle-vibration coupling, and quasi-particle residue at the pole, typical 
of NFT, allows theory to make predictions which are more accurate than the 
basic parameters entering the calculation, as a result of a delicate process of 
cancellation of errors. 

4 Conclusions 

Summing up, we have found that the small difference existing between the 
matrix elements of the bare viow-k interaction and those of the Gogny in- 
teraction leads to important differences between the pairing gaps associated 
with the two forces {Acogny ~ '^^low-k, for rrik ~ 0.7m). This difference is 
removed and eventually overwhelmed by including the attractive contribu- 
tion coming from the induced interaction Vind arising from the exchange of 
surface vibrations, which acts only on a rather small energy range around ep- 
We note that the pairing gap in ^^''Sn is not much changed, including also the 
effect of spin fluctuations |22I; which instead give the dominant (repulsive) 
contribution to the induced interaction in neutron matter. 

The 25% excess displayed by the resulting total pairing gap with respect 
to the experimental value is corrected by considering the corresponding quasi- 
particle strength at the Fermi energy. An important question which remains 
to be addressed quantitatively is the ability viow-k + Vind as well as Vcogny 
have to describe in detail the isotopic and nuclear structure dependences 
displayed by Aexp- This subject will be addressed in a future publication. 

Within this context, we note that the dependence of the pairing gaps 
generated by vcogny and by viow-k + Vind on temperature and rotational fre- 
quency is expected to be quite different, as vcogny and viow-k are essentially 
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independent of these parameters, while Vind is quite sensitive to them. 

It is also worth mentioning that a recent study in symmetric nuclear 
matter revealed a very strong extra attraction coming from the induced in- 
teraction [221 ■ Nuclear matter and finite nuclei may, however, show only 
qualitative similarities in this case and a quantitative link may be difficult 
to establish. 



We thank L. Coraggio for providing us with the matrix elements of the 
Viow-k interaction. 
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Appendix A 



Simple estimate of the pairing gap 



In what follows we consider some of the consequences the particle- vibration 
coupling has on the pairing correlations of particles moving in a single j-shell 
interacting through a bare nucleon-nucleon pairing potential with constant 
matrix elements G 

For this simple model, the value of the occupation numbers Ui, and 
must be the same for all the 2j + 1 orbitals. In particular, the occupation 
probability for the case when the system is occupied with N = Q particles 
(half filled shell), where 

(A.1) 

is 




U = \ll-i^ = \ll- (A.3) 
Consequently, the pairing gap is given by the relation 

A = GJ2U.V^ = ^- (A.4) 

Because the density of levels is proportional to the a;— mass (m^^), the effective 
(dressed) degeneracy can be written 

^e// = ^, (A.5) 

in terms of the quasiparticle strength at the Fermi energy = {m^/m)~^ 
where m^^ = m(l + A), A being the mass enhancement factor. In the case of 
nuclei this dimensionless quantity (which measures the strength with which 
nucleons couple to low-lying collective vibrations) is of the order of 0.5, closer 
to the strong than to the weak coupling situation of BCS (in which case 
A ~ A) i2I|. 



13 



Because of their coupling to vibrations, nucleons spend part of the time in 
more comphcated configurations than pure single-particle states. The quan- 
tity (quasiparticle strength |9j) = (1 -|- A)~^ ~ 0.7 measures the content 
of single-particle strength present in levels around the Fermi energy avail- 
able for nucleons to interact through a two-body interaction, in particular 
through through a pairing force. Consequently, due to the self-energy pro- 
cesses arising from the particle-vibration coupling phenomenon, the pairing 
strength becomes GZj'. 

The exchange of vibrations between pairs of nucleons moving in time 
reversal states close to the Fermi energy gives rise to an effective pairing 
interaction of strength g^^yZj' ^ where g^^^ stands for the induced pairing 
interaction controlled by the particle-vibration coupling vertices. Taking 
these effects into account, one can write 

A = |^(G + (7p_,)r] (A.6) 

Identifying G with viow-k and g^^^ with vind-, one finds from Fig. IHl that 

^(G + (?p_,)fi^2MeV. (A.7) 

As seen from Eq. ()A.6|) . the number to be compared with experimental find- 
ings is 

A^en = Z^ X 2MeV ^ 1.4MeV. (A.8) 

Note that in the above discussion we have not considered the errors to be ex- 
pected in estimating both 5fp_^ and Z^. This subject is taken up in Appendix 
B. 
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Appendix B 



Simple estimate of A and 



From Fig. [T] it is seen that the average values of Vind and viow-k at the 
Fermi energy are < Vind > = -0.15 MeV and < viow-k > = -0.16 MeV. Because 
Vind is of the order of J2l Pl/^^-^l, the matrix elements of Vind will display 
an error equal to twice the average error displayed hj Pl plus that displayed 
by DENl- Assuming this error to be 20% for both quantities (a rather 
extreme negative situation) one concludes that < Vmd >= 0.15±0.09 MeV. 
Consequently < Vtot >=< Vind > + < viow-k >= — 0.31±0.09 MeV. Because 
Atotai = 2 MeV for < Vtot > = -0.31 MeV, one obtains At^tai = (2 ± 0.6) 
MeV. 

We now proceed to the calculation of A = A^(0) < Vind >, where A^(0) is 
the density of single-particle levels around the Fermi energy for one spin ori- 
entation and one type of nucleon. Making use of the value iV(0) = 3.4 MeV~^ 
appropriate for ^^°Sn, one gets A = 0.5 ± 0.3. Consequently A^en = Z^^Atot 
= 1.35 ± 0.15 MeV, where Z^; = (1 -|- A)^^. We thus note a conspicuously 
smaller error for A^en than for Vmd- This is because the larger is A (induced 
pairing contribution), the smaller is the content of single-particle strength 
and thus of Z^, associated with levels lying close to the Fermi energy, and 
viceversa. 
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Figure 1: The nucleus ^^'^Sn. Diagonal pairing matrix elements of the induced 
interaction (upper panel, solid diamonds) and of the Gogny force (lower 
panel, solid circles), displayed as a function of the single-particle energy 
of the state z/ calculated using the bare nucleon mass and the single-particle 
wavefunctions of a Woods-Saxon potential with standard parameters (depth 
Vo= -49 MeV, diffusivity a =0.65 fm, and radius Rq = 6.16 fm), including the 
spin-orbit term, parametrized according to ref. [201 ■ Also shown by means 
of vertical lines is the position of the Fermi energy, ep = —9.1 MeV. Note 
the different scale in the two figures. 
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Figure 2: Schematic model for the nucleus ^^°Sn. The semiclassical pairing 
matrix elements of the Gogny force (solid line) are compared to the quan- 
tal matrix elements for the case in which the single-particle wavefunctions 
have been calculated making use of a harmonic oscillator potential without 
the spin-orbit term (and = m). The diagonal and non-diagonal matrix 
elements are denoted by filled and open circles respectively. The crosses 
correspond to the weighted averages of the matrix elements within a shell. 



17 




-40 -30 -20 -10 

[MeV] 



Figure 3: The nucleus ^^°Sn. The same as Fig. |21 (i.e. also with = m) 
but for the fact that the single-particle wavefunctions have been calculated 
making use of the Woods-Saxon potential including the spin-orbit term al- 
ready used for Fig. Q where the diagonal matrix elements have already been 
shown. The non-diagonal matrix elements are plotted here at an energy 
which is the average between the energies of the initial and final states. 
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Figure 4: The matrix elements of the Gogny (soUd curve) and of the Viow-k 
(dashed curve) interactions are plotted as a function of momentum. 
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Figure 5: The nucleus ^^°Sn. The semiclassical matrix elements of the 
induced interaction, calculated according to Eq. (dash-dotted curve), are 
compared with the matrix elements of the Gogny force (solid curve, cf. Fig.EI) 
and with those of the viow-k interaction (dashed curve). Calculations are 
performed with = m and with the same Woods-Saxon potential used in 
Figs. [Hand El 
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Figure 6: The nucleus ^^°Sn. The semiclassical matrix elements of the in- 
duced interaction, as a function of the energy of the single-particle levels 
(calculated with = m) associated with the multipolarities L = 2, 3, and 
4 of the formfactor (see Eq.(PQ)), are shown by a dashed, dotted and solid line 
respectively. The sum of the three contributions is displayed by means of a 
dash-dotted line (cf. Fig. EJ. 
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Figure 7: Pairing gaps calculated in neutron matter as a function of the Fermi 
momentum, obtained with the Gogny interaction (solid line), the Argonne 
vu potential (dash-dotted line) and the viow-k potential (dashed line). The 
bare effective mass has been used in the calculation. 
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Figure 8: State-dependent pairing gaps of ^^°Sn calculated with a Woods- 
Saxon potential (with depth Vq = -64 MeV, diffusivity a = 0.65 fm and radius 
i?o = 6.17 fm), as a function of the single-particle energy. The /c— mass 
was set equal to 0.7m. The Fermi energy is ep —- 8.6 MeV. Sohd triangles 
(open squares) display the results of a HFB calculation with the Gogny in- 
teraction, with quantal (scmiclassical) matrix elements. The solid diamonds 
refer instead to a HFB calculation using the semiclassical matrix elements of 
the viow-k potential. 
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Figure 9: State-dependent pairing gap of Sn obtained making use of the 
Woods-Saxon potential used in calculating the results displayed in Fig. |Hl(i.e. 
with ruk = 0.7m) and different interactions. Diamonds and triangles show 
again the gaps obtained with the viow-k and Vcogny interactions, respectively. 
Stars show the gap associated with the matrix elements obtained summing 
the matrix elements of viow-k and Vmd- The latter have been evaluated using 
Eqs. ® and ^ 
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